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Abstract. We undertake a detailed numerical study of the twin phenomenon of stochastic and vibrational 
resonance in a discrete model system in the presence of bichromatic input signal. A two parameter cubic 
map is used as the model that combines the features of both bistable and threshold settings. Our analy- 
sis brings out several interesting results, such as, the existence of a cross over behavior from vibrational 
to stochasic resonance and the possibility of using stochastic resonance as a filter for the selective detec- 
tion/transmission of the component frequencies in a composite signal. The study also reveals a fundamental 
difference between the bistable and threshold mechanisms with respect to amplification of a multi signal 
input. 

PACS. 05.10-a Computational methods in statistical and nonlinear physics - 05.40 Noise 

1 Introduction noisy systems has increased due to the applications in the 

modelling of a great variety of phenomena of relevance in 

During the last two decades, investigation of signal pro- physics, chemistry and lifesciences[l-3]. 

cessing in nonlinear systems in the presence of noise has T on ii r r ± x. i • 

° J r In bK, the response ol a nonlinear system to a weak m- 

revealed several interesting phenomena, the most impor- , . n . . . c ,-, . , .,, , 

° r ' r put signal is significantly increased witn appropriate tun- 

tant being stochastic resonance f SR.) [11 . It can roughly be <■ • , rn \m 1,1 u u • i 

to ° J mg oi the noise intensity [lj. When a subthreshold signal 

considered as the optimisation of certain system perfor- r /,\ . • , , ,. . . ,i 

* 1 [tj is input to a nonlinear system g together witn a noise 



mance by noise. The interest in the study of nonlinear 



C(i), if the filtered output 0(t) = g(I(t) + £(*)) shows 
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of J(t), then SR is said to be realised in the system. The If C t h is the threshold at which deterministic switching 

mechanism first used by Benzi, Nicolis etc [4,5] to explain (with noise amplitude E = 0) is possible, then well to well 

natural phenomena is now being used for a large variety of switching due to SR occurs when 
interesting applications like modelling biological and eco- 

- C th > (ZStriiut + E((t)) > C th (3) 

logical systems [6], lossless communication purposes etc 

[7]. Apart from these, it has opened up a vista of many that is, twice in one period of the signal T = 2tt/w. 

related resonance phenomena [8] which are equally chal- The characterisation of SR in this case is most com- 

lcnging from the point of view of intense research. The monly done by computing the signal to noise ratio (SNR) 

most important among these, which closely resembles SR, from the power spectrum of the output as 



SNR = mog w {S/N)dB (4) 



is the vibrational resonance (VR) [9], where a high fre- 
quency forcing plays the role of noise and amplify the re- 

sponse to a low frequency signal in bistable systems. In where N is the average background noise around the signal 



VR, analogous to SR, the system response shows a bell S - If SR occurs in thc s y stcm > thcn thc SNR § ocs through 

shaped resonant form as a function of thc amplitude of a maximum giving a bell shaped curve as E is tuned, 
thc high frequency signal. In this work, we try to capture SR has also been observed in systems with a single 

numerically some interesting and novel aspects of SR and stable state with an esca P e scenerio. These threshold sys- 

VR, using a simple discrete model, namely, a two param- tems > in the simplest case, can be modelled by a piecewise 

eter bimodal cubic map. lincar svstem or ste P function 



g(u) = -1 u < C th 



In thc early stages of the development of SR, most 
of the studies were done using a dynamical setup with 
bistability, modelled by a double well potential. Here SR 
is realised due to the shuttling between the two stable 

states at the frequency of the subthreshold signal with = +1 u > Cth 

the help of noise. Thus if the potential is 



V(X) = -aX 2 /2 + bX 4 /4 (1) 



The escape with the help of noise is followed by reset- 
ting to the monostable state. In this case, a quantitative 

. £ ■ i j • ,i j • , characterisation is possible directly from the output, but 

m thc presence of a signal and noise, the dynamics can be ^ J ^ ' 

, , .,, , only in terms of probabilities. If t n are the escape times, 

modelled by an overdamped oscillator J * r ■> 

the inter spike interval (ISI) is defined as T n — t n+ \ — t n 
X = aX — bX 3 + ZSinut + E((t) (2) and m(T n ) is thc number of times the same T„ occurs. 
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pled noisy oscillators [11] and two coupled anharmonic os- 
cillators [12] in the bistable setup. In the threshold setup, 
VR has been shown to occur only in one system, namely, in 
the numerical simulation of the FitzHugh-Nagumo (FHN) 
model in the excitable regime along with the experimen- 
tal confirmation using an electronic circuit [13]. Here we 
show for the first time the occurance of VR in a discrete 
system, both in bistable and threshold setups. 

There are situations where SR and VR are to be opti- 
mised by adapting to or designing the dynamics of the sys- 
tem. This is especially relevant in natural systems where 
the noise is mostly from the environmental background 
and therefore not viable to fine tuning. Similarly, depend- 
ing on the context or application, the nature of the signal 
can also be different, such as, periodic, aperiodic, digital, 
composite etc. The classical SR deals with the detection of 
a single subthreshold signal immersed in noise. However, 
in many practical situations, a composite signal consist- 
ing of two harmonic components in the presence of back- 
ground noise is encountered, as for example, in biologi- 
cal systems for the study of planktons and human visual 
cortex [14], in laser physics [15] and in accoustics [16]. 
Studies involving such bichromatic signals are also rele- 
vant in communication, since one can address the question 
of the carrier signal itself amplifying the modulating sig- 
nal. Moreover, two frequency signals are commonly used 
in multichannel optical communication systems based on 
wave length division multiplexing (WDM) [17]. But only 
very few studies of SR have been carried out using the 
bichromatic signals to date [18-21], each of them pertain- 



Fig. 1. The bifurcation structure of the cubic map given in 
eq.(5) corresponding to four values of the parameter a, to show 
the bistability in periodic and chaotic attractors with the basin 
boundary at X = 0. It is clear that the bistability just begins 
at a — 1.0 and ends at a = 2.6. Note that the range of b values 
for bistability is different for periodic and chaotic attractors. 

For SR to be realised in the system, the probability p n = 
m(T n )/N (N is the total number of escapes) has to be 
maximum corresponding to the signal period T at an op- 
timum noise amplitude. 

In the case of VR, the high frequency forcing takes 
the role of noise to boost the subthreshold signal. The 
system is then under the action of a two frequency sig- 
nal, one of low frequency and the other of high frequency, 
with or without the presence of noise. Most of the studies 
in VR also have been done using the standard model of 
overdamped bistable oscillators [10]. But recently, VR has 
been shown to occur in a spatially extended system of cou- 
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ing to some specific dynamical setups and with continuous 
systems. This is one of the motivations for us to under- 
take a detailed numerical analysis of SR and VR with such 
signals using a discrete model. It is a discretised version 
of the overdamped bistable oscillator, but with the added 
feature of an inherent escape mode also. Hence it can func- 
tion in both setups, bistable and threshold, as a stochastic 
resonator. 

There are many situations, especially in the biologi- 
cal context, such as, host-parasite model, virus-immune 
model etc., where discrete systems model the time devel- 
opments directly. They can behave differently, especially 
in the presence of high frequency modulation and back- 
ground noise. The benefit of high frequency forcing has 
been studied in the response of several biological phenom- 
ena^]. High frequency stimulation treatments in Parkin- 
son's disease and other disorders in neuronal activity have 
also been reported[23]. Moreover, it is also known that op- 
timum high frequency modulation improves processing of 
low frequency signal even in systems without bistability 
where noise can induce the required structure[ll]. Thus a 
study of resonance phenomena in discrete systems can lead 
to qualitatively different results having potential practical 
applications. We do observe some novel features which 
have not been reported so far for continuous systems. 



■1.5 



Fig. 2. Parameter space plot in (a, b) plane for the cubic 
map given by eq.(5). Regions of periodic attractors with pe- 
riod l(p=l), period 2(p=2) etc. and that of chaotic attractors 
and escape are shown in different colours. Region within the 
quadrilaterals (a > 1.0) corresponding to p = l,p = 2 and 
chaos represent the respective bistable regions. 

SR and VR with the cubic map as a threshold system. 
Results and discussions are given given in §6. 

2 The model system 

The model used for our analysis is a two parameter cubic 
map given by 



X n+1 = f(X n )=b + aX n -Xl 



(5) 



Our paper is organised as follows: In §2, the model sys- It is a discrete version of the Duffing's double well po- 

tem used for the analysis is introduced. In §3, we study SR tential and is the simplest possible nonlinear discrete sys- 

in the model numerically with bichromatic signal treating tern that combines the desirable features of bistable and 

it as a bistable system. Numerical and analytical results of threshold setups. Hence it is extremely useful in the study 



VR in the bistable setup are presented in §4. §5 discusses of resonance phenomena caused by noise or high frequency 
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signal. Similar systems in the continuous cases include the where we choose ((n) to be a gaussian noise time series 

FHN model[23] for neuronal firing with two widely differ- with zero mean and F(n) is the periodic signal sampled 

cnt time scales. in unit time step. The amplitude of the noise and the 

The system has been studied in great detail both ana- signal can be varied by tuning E and Z respectively. It 

lytically and numerically and has been shown to possess a can be shown that in the regime of chaotic bistable at- 

rich variety of dynamical properties including bistability tractors (a = 2.4,6 = 0.01), a subthreshold input sig- 

[24]. In particular, if a\ is the value of the parameter at nal can be detected using the inherent chaos in the sys- 

which f'(Xi, a\,b) = 1, then for a > ai, there is a window tern without any external noise(E=0). Taking the signal 

in 6, where bistability is observed. The bistable attractors F(n) = ZSin{2irvn), with Z = 0.16, the system shows 

are clearly separated with X > being the basin of one SR type behavior for an optimum range of frequencies 

and X < that of the other. It is found that for a in the as shown in Fig. 3, where the output SNR is plotted as a 

range 1.0 < a < 2.6, the system has periodic states and function of the frequency v. It implies that a subthresh- 

chaotic states. As a increases from 1.0, the periodicity of old signal can be detected by passing through a bistable 

the bistable attractors keep on doubling while the width system making use of the inherent chaos in it without the 

of the window around b decreases correspondingly. For help of any external noise. This phenomenon is known as 

a = 2.4, two chaotic attractors co-exist in a very narrow deterministic resonance. In the regime of periodic bistable 

window around b. For a < 1.0, the system has a monos- attractors (a = 1.4,6 = 0.01) with a single subthreshold 

table period 1 state and for a > 2.6, there is merging of the signal, the system shows conventional SR as well as chaotic 

chaotic states followed by escape. All these can be clearly resonance (CR), and using this model, we have recently 

seen from Fig.l where bifurcation structure of the cubic reported some new results including enhancement of SNR 

map corresponding to four a values are shown. A detailed via coupling [25]. 
stability analysis fixes the different asymptotic behavior 

in its parameter plane (a, b) as shown in Fig.2. Regions of 3 SR with trichromatic input signal 

periodicities (p = l,p = 2) etc., chaos and escape in the 
parameter plane can be clearly seen. The quadrilateral re- 
gions in the area marked p = l,p= 2 and chaos represent 
bistability in the respective attractors. 

The system when driven by a gaussian noise and a 
periodic signal becomes 



In this section, we undertake a numerical study of SR 
using the cubic map in the bistable regime. The input 
signal F(n) is taken as a superposition of at least two 
fundamental frequencies, v\ and v 2 - Before going into the 
numerical results, we briefly mention some already known 
results for multisignal inputs in the case of over damped 
X n+1 = b + aX n -X 3 n + E((n) + ZF(n) (6) bistable potential. 
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Fig. 3. Variation of SNR with frequency for the bistable cubic 
map in the chaotic regime (a = 2.4) with the noise amplitude 
E = 0. Note that a subthreshold signal = 0.16) can be 
detected for an intermediate range of frequencies, without any 
external noise. 

The most popular theory for the analytical description 
of SR is the linear response theory(LRT) [26-29]. Accord- 
ing to this theory the response of a nonlinear stochastic 
system X(t) to a weak external force F(t) in the asymp- 
totic limit of large times is determined by the integral 
relation [26] 



R(t - T)F(r)d7 

-OO 



(7) 



where < Xq > is the mean value of the state variable 
for F(t) — 0. Without lack of generality, one can set 
< X >= 0. The function R(t) in (7) is called the re- 
sponse function. For a harmonic signal, the system re- 
sponse can be expressed through the function R(u>) which 



is the Fourier transform of the response function: 



X(t) = Z\R{uj)\Sin{2irvt + tp) 



(8) 



where ip is a phase shift. 

The LRT can be naturally extended to the case of mul- 
tifrequency signals. Let the signal F(t) be a composite 
signal of the form: 



F(t) = Z^Sin{2irv k t) 



(9) 



fc=i 



where v k are the frequencies of the discrete spectral com- 
ponents with the same amplitude Z. Then according to 
LRT, the system response can be shown to be [29] 

n 

X{t) = zY J \R{^k)\Sin{2-Kv k t + ^ k ) (10) 
fc=i 

which contains the same spectral components at the in- 
put. We now show numerically that the same is true for a 
discrete system as well, in the bistable setup. 

For the remaining part of the paper, we fix a — 1.4 
and b = 0.01 in the region of bistable periodic attractors 
and ((n) is a gaussian time series whose amplitude can be 
tuned by changing the value of E. The system is driven 
by a composite signal consisting of a combination of two 
frequencies v\ v%, given by 



ZF(n) = Z{Sin2iruin + Sin2~KVin) 



(11) 



The value of Z is fixed at 0.16 so that the signal is well 
below threshold. It should be noted that because of the 
iteration with unit time step (ze, n = 1, 2, 3, ...), the avail- 
able range of frequencies is limited to Vk = (0,0.5). We 
use different combination of frequencies (1/1,1/2) for nu- 
merical simulation in the presence of noise by tuning the 
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Fig. 4. Power spectrum for the bistable system (6) 
with a composite signal consisting of (a) 2 frequencies 
(yi,V2) = (0.125,0.05) and (6) 3 frequencies (vi, v%, v$) — 
(0.125,0.075,0.2), with Z = 0.16 and E = 0.4 in both cases. 
The system parameters are a = 1.4 and b — 0.01. 

value of E. For convenience, v\ is fixed at 0.125 and is 
varied from 0.02 to 0.5 in steps of 0.01. For each selected 
combination, the output power spectrum is calculated us- 
ing the FFT algorithm for different values of E. A typical 
power spectrum for (^1,^2) = (0.125,0.05) is shown in 
Fig. 4(a). The procedure is repeated with F(n) consisting 
of a combination of 3 frequencies and a power spectrum 
for (^,^2,^3) = (0.125,0.075,0.2) is shown in Fig.4(b). 
To compute the power spectrum, only the inter-well tran- 
sitions are taken into account and all the intra-well fluctu- 
ations are suppressed with a two state filtering. It is clear 
that, only the fundamental frequencies present in the in- 
put are enhanced. 
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Fig. 5. Variation of SNR with noise amplitude for the 2 fre- 
quencies shown in Fig. 4(a), with the filled circle connected by 
solid line representing the higher frequency v\. 



We now concentrate on a combination of 2 frequen- 
cies and compute the two important quantifiers of SR, 
namely, the SNR and the Residence Time Distribution 
Function(RTDF). For the frequencies in Fig.4(a), the SNR 
is computed from the power spectrum using equation (4) 
for a range of values of E and the results are shown in 
Fig. 5. The RTDF measures the probability distribution of 
the average times the system resides in one basin, as a 
function of different periods. If T is the period of the ap- 
plied signal, the distribution will have peaks corresponding 
to times (2n + l)T/2, n=0,l,2...For the system (6) with 
(^1,^2) = (0.125,0.05), the results are shown in Fig. 6. 
Note that there are only peaks corresponding to the half 
integer periods of the two applied frequencies. 
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Fig. 6. The RTDF for the system (6), with a composite signal 
consisting of 2 frequencies v\ — 0.125 and vi — 0.05 and am- 
plitude Z — 0.16. The noise amplitude is put at the optimum 
value 0.4. Note that the peaks are synchronised with half the 
periods corresponding to the input frequencies. 

The above computations are repeated taking various 
combinations of frequencies (^1,^2), both commensurate 
and non-commensurate. For a fixed combination of (y\ , ^2) , 
the calculations are also done by changing the signal am- 
plitude Z of one and both signals. Always the results 
remain qualitatively the same and only the fundamen- 
tal frequencies present in the input are amplified at the 
output. If Z becomes very small (< 0.1) compared to 
the noise level, then the phenomenon of SR disappears 
altogether and the signal remains undetected in the back- 
ground noise. 

In all the above computations, we used additive noise, 
where the noise has been added to the system externally. 



But in many natural systems, noise enters through an in- 
teraction of the system with the surroundings, that is, 
through a parameter modulation, rather than a simple 
addition. Such a multiplicative noise occurs in a variety 
of physical phenomena [30] and can, in principle, show 
qualitatively different behavior in the presence of a peri- 
odic field [31]. To study its effect on the bistable system, 
equation (6) is modified as 

X n+1 = b + a(l+ E((n))X n - X* + ZF(n) (12) 

The noise is added through the parameter a which deter- 
mines the nature of the bistable attractors. With a = 1.4 
and b = 0.01, the system is now driven by a signal of single 
frequency v\ = 0.125 and a multisignal with 2 frequencies 
(^1, ^2), with Z = 0.16. The power spectrum for single fre- 
quency and multiple frequencies are shown in Fig. 7(a) and 
(b) respectively. The corresponding SNR variation with 
noise E are shown in Fig. 8 (a) and (b). Note that the re- 
sults are qualitatively identical to that of additive noise, 
but the optimum SNR and the corresponding noise ampli- 
tude are comparitively much higher in this case. Thus our 
numerical results indicate that a bistable system responds 
only to the fundamental frequencies in a composite signal 
and not to any mixed modes. 

4 VR: Numerical and analytical results 

In order to study VR in the system, the input periodic 
driving is modified as 

ZF{n) = Z{l)Sin{2-KV X n) + Z{2)Sin{2Trv 2 n) (13) 
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Fig. 7. The power spectrum for the bistable system (12), with 
the signal F(n) consisting of (a) one frequency v\ = 0.125 
and (6) 2 frequencies v\ — 0.125 and V2 — 0.3. The parameter 
values used are Z — 0.16, E = 1.0, a = 1.4 and b — 0.01. 

where vi is the low frequency signal which is fixed at 0.125 
with its amplitude Z(l) at the subthreshold level 0.16. It 
is added with a signal of high frequency v%(> ^i) whose 
amplitude Z(2) is tuned to get VR. We have used a num- 
ber of values for v-i over a wide range from 0.15 to 0.5 to 
study VR in the system. A small amount of noise is also 
added as in eq.(6) which represents tiny random fluctua- 
tions present in all practical systems. 

By taking — 0.4 and E = 0.12, the system is iter- 
ated by tuning the high frequency signal amplitude. The 
variation of SNR for the signal computed from the output 
time series as a function of Z(2), clearly indicates VR in 
the system. The procedure is repeated by changing E to 
study the influence of noise level on VR. The results are 
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Fig. 8. Variation of SNR with noise for the system(12), for (o) 
single frequency v\ = 0.125 and (b) 2 frequencies v\ — 0.125 
(filled circles) and = 0.3. Note that the optimum SNR of V\ 
increases by about 5 dB when a second signal V2 is added. 

shown in Fig. 9 for 3 values of E. It is found that as E de- 
creases, the optimum SNR shifts towards the higher value 
of Z(2). This result has been reported earlier in the case 
of continuous systems also[10,12]. 

The computations are repeated by taking different val- 
ues of V2 in the range 0.15 to 0.5 by fixing E — 0.12. The 
SNR as a function of Z(2) for four values of Vi are shown 
in Fig. 10. Two results are evident from the figure. The 
optimum SNR is independant of the high frequency when 
{i>2 — v\) is small. But for v% >> v\, the optimum SNR 
decreases appreciably. This implies that the effectiveness 
of the high frequency in producing VR is reduced as it be- 
comes too large compared to the low frequency signal. It 
can also be seen that the optimum value of Z(2) depends 
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Fig. 9. Variation of SNR corresponding to the frequency 
v\ — 0.125 as a function of the amplitude Z(2) of the high 
frequency signal with = 0.4, showing VR in the bistable cu- 
bic map with a bichromatic input signal given by eq.(13). The 
three curves correspond to different values of noise amplitude, 
namely, E = 0.06 (circle connected by solid line), E — 0.12 
(triangle connected by dotted line) and E — 0.16 (squares con- 
neted by dashed line). The optimum values of Z(2) shifts to- 
wards the left as E increases. 

on v-i. Infact, it increases with vi. Thus, a higher value of 
vi requires a larger amplitude and becomes less effective 
as i*2 >> v\. 

One reason for this is that, due to the unit time sam- 
pling of the signals, the fluctuation in the signal values 
will be more as the frequency increases. Hence the av- 
erage will be closer to zero, higher the frequency, which 
makes it necessary to have larger amplitude. But once the 
optimum amplitude becomes too large, escape becomes 



Fig. 10. Variation of SNR for v\ as a function of Z(2) cor- 
responding to four values of vi which are 0.18 (filled circle 
connected by solid line), 0.32 (triangles connected by dotted 
line), 0.40 (squares connected by dashed line) and 0.46 (open 
circles connected by solid line). As ^2 >> v\, it becomes less 
effective in producing VR as indicated by the reduction in the 
optimum value of SNR. 

possible in all time scales making it less effective. It can 
also be shown that the high frequency forcing can change 
the effective value of the parameter a so that the map can 
come out of the bistable window. This is discussed in more 
detail below. 

The dependance of optimum Z(2) on v-i also implies 
that the system shows the so called bona fide resonance[32,33] . 
This is shown in Fig.ll, where the SNR is plotted as a 
function of the high frequency for two values (0.5&0.7) of 
Z(2), with E fixed at 0.12. 
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We now briefly discuss some analytical results for VR Thus, for a given choice of a, there is an upper limit for 

in the cubic map and show that the high frequency forcing (J^) for retaining the bistability in the system. The actual 

can change the effective value of the parameter a. Here we values may get modified in the presence of noise. Also, this 

consider the case v 2 » v\, in the noise free limit, E = 0. mechanism provides two parameters B and w 2 that can 

In parallel with the theory developed to explain VR in be tuned in a mutually compramisable manner to obtain 

continuous systems [34], we analyse the effect of the widely VR or even SR with added noise. Moreover, even when 

differing frequencies for the system under study: a > 2.6, a* can be < 2.6 and hence this increases the vir- 

~ . tual window of bistability in the system providing greater 

X n +i =b + aX n - Xl + ASmum + BSmu 2 n (14) J J 66 



range for applicability. 



where Wi ; 2 = 2-kv\^- Taking A = 0, we look for a solu- 
tion^] 



BSinu 2 n 5 SR and VR in the threshold setup 

X n =Y n 2 (15) 

While the first term Y n varies significantly only over time As said earlier, the domains of the bistable attractors in 

of the order n, the second term varies rapidly within an the cubic map are clearly separated with the boundary 

iteration and hence can be averaged. Putting (15) in (14) X = 0. Hence the cubic map can also be considered as a 

and averaging, we get: nondynamical threshold system with a single stable state 

3Y n B 2 having a potential barrier. Here the system generates an 



Y n+1 = b + aY n - Yl - (16) 



-2 



That is, 



where 



Y n+1 = b + a*Y n - V„ 3 (17) 



output spike only when the combined effort of the signal 
and the noise pushes it across the potential barrier (atX = 
0) in one direction: 



3B 2 [E((n) + ZF{n)\ > C th (21) 

fl * = a -^f (18) 

Tn, 4-u a v c 4-u u- u r r • • *- j- ^ It is then externally reinjected back into the basin. The 

1 hus the effect of the high frequency forcing is to reset the J J 

parameter a as a*. Hence only for the choice of B and lu 2 out P ut thus consists of a series of s P ikes similar to a ran " 

that retains a* in the bistable window (1.0 < a* < 2.6), do dom telegraph process. The study of SR and VR in such 

• uj-ii- u i, 4- j--u i -c • 1 systems assumes importance in the context of biological 

we expect shuttling behavior at the low frequency signal ■' v ° 

c a —L a at 4-t, 1 r -4- r u applications and in particular the integrate and fire mod- 
tor A ^ 0. Now the lower limit for a becomes rv r t> 

oto2 els of neurons where SR and VR have been firmly estab- 
o>1.0+— — (19) 

2tJ 2 lished[3,35]. 

which gives The computations are done using equation (6), ini- 
B 2 

~3 ( ' a ~ 1 ' ' ) ^ tially with a single frequency signal, F(n) = Sin{2'KV 1 n). 
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Fig. 11. Evidence for bona-fide resonance in the cubic map. 
Variation of SNR corresponding to v\ = 0.125 as function of 
the high frequency ^2 for two values of its amplitude, namely, 
Z{2) = 0.5 (circles) and Z(2) = 0.7 (squares). The amplitude 
of v\ is fixed at 0.16 and E = 0.12. 

We start from an initial condition in the negative basin 
and when the output crosses the threshold, X = 0, it is 
reinjected back into the basin by resetting the initial con- 
ditions. This is repeated for a sufficiently large number 
of escapes and the ISIs are calculated. The ISIs are then 
normalised in terms of the periods T n and the probability 
of escape corresponding to each T n is calculated as the 
ratio of the number of times T n occurs to the total num- 
ber of escapes. The whole procedure is repeated tuning 
the noise amplitude E. It is found that the ISI is synchro- 
nised with the period of the forcing signal for an optimum 
noise amplitude (Fig. 12), indicating SR for the frequency 
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Fig. 12. The probability of escape for the system (6) as a func- 
tion of different periods, when it is used as a threshold system 
with the signal F(n) having a single frequency v\ — 0.125 and 
a subthreshold amplitude Z = 0.16. The escape probability 
is plotted corresponding to four different noise amplitudes E. 
Note that the probability corresponding to the signal period 
(T = 8) passes through an optimum value for E — 0.5. 

The calculations are now repeated by adding a second 
signal of frequency v<i and amplitude same as that of v\ . 
Again different values of v 2 in the range 0.02 to 0.5 are 
used for the calculation. It is then found that apart from 
the input frequencies v\ and 1^2, a third frequency, which 
is a mixed mode is also enhanced at the output, at a lesser 
value of the noise amplitude. To make it clear, the ampli- 
tude Z of both v\ and v 2 are reduced from 0.16 to 0.08, 
so that they become too weak to get amplified. The re- 
sults of computations are shown in Fig. 13 and Fig. 14, for 
a combination of input signals (^1,^2) = (0.125,0.033). 
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Fig. 13 represents the probability of escape corresponding 
to different periods, for the optimum value of noise. It 
is clear that only a very narrow band of frequencies dv 
around a third frequency « 0.045 (corresponding to the 
period T ~ 22s), are amplified at the output. Note that 
this frequency is absent in the input and corresponds to 
(y\ — vi)j1. This is in sharp contrast to the earlier case 
of a bistable system. The variation of escape probability 
corresponding to this frequency as a function of noise am- 
plitude is shown in Fig. 14. 

This result can be understood as follows: When two 
signals of frequencies v\ and v-i and equal amplitude Z 
are superposed, the resulting signal consists of peaks of 
amplitude 2Z repeating with a frequency (y\ — v-i)j2 in 
accordance with the linear superposition principle: 

Sin{2nvit) + Sin{2nv2t) = 2Sin{2nvj r t)Cos{2 , nv^t) 

(22) 

where v+ = {vi+V2)/2 and V- = (v\—V2)/2. For athresh- 
old system, the probability of escape depends only on the 
amplitude of the signal which is maximum corresponding 
to the frequency V- . But in the case of a bistable system, 
the signal is enhanced only if there is a regular shuttling 
between the wells at the corresponding frequency. This 
is rather difficult for the frequency V- because, its am- 
plitude is modulated by a higher frequency z/+. This re- 
sult has been checked by using different combinations of 
frequencies (z-'i,^) and also with different amplitudes. It 
should be mentioned here that this result reveals a fun- 
damental difference between the two mechanisms of SR 
and is independant of the model considered here. We have 
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Fig. 13. The escape probability corresponding to the optimum 
noise amplitude with the signal F(n) comprising of 2 frequen- 
cies v\ = 0.125 and v-2 = 0.033. The amplitudes of individual 
signals are put much below the threshold value required for SR. 
Note that only the difference frequency {y\ — i/2)/2 is amplified. 

obtained identical results with a fundamentally different 
model showing SR, namely, a model for Josephson junc- 
tion and has been discussed elsewhere [21]. 

To study VR in the threshold set up, we once again 
consider the input signal as in eq.(13) with the low fre- 
quency v x fixed at 0.125 and Z(l) = 0.16. With E = 0.12 
and V2 = 0.32, the computations are performed as be- 
fore by tuning the high frequency amplitude Z(2). The 
results are shown in Fig. 15, which clearly indicates VR in 
the threshold set up. The computations are repeated for 
different values of V2 and here also it is found that the 
optimum escape probability becomes less for > 0.45 as 
in the bistable case. 
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Fig. 14. Variation of the probability of escape correspond- 
ing to the frequency shown in Fig. 13, as a function of noise 
amplitude. 



It is interesting to note that when the high frequency 
forcing and noise are present simultaneously, only one of 
them dominates to enhance the signal. In the regime of 
VR, the noise level has to be small enough and vice versa 
for SR. When both are high, the transitions become ran- 
dom and the signal is lost. We find that the two regimes 
of (VR and SR) can be distinguished in terms of an initial 
time delay for the system to respond to the high frequency 
or noise, called the response time. The response time (r) 
decreases sharply in the regime of SR (ie, when signal is 
boosted with the help of noise), thus indicating a cross 
over behavior between the two regimes. This is shown in 
Fig. 16, where r is plotted as a function of E for the thresh- 
old setup. The computations are done in such a way that 
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Fig. 15. Escape probability as a function of different periods 
corresponding to four values of the signal amplitude of the high 
frequency (1/2 = 0.32), showing VR for the low frequency signal 
v\ for the cubic map in the threshold setup. Other parameters 
are a = 1.4, b = 0.01, E = 0.12 and Z(l) = 0.16. 

as E is increased, Z(2) is decreased correspondingly to get 
optimum response. 

6 Results and discussion 

In this paper we undertake a detailed numerical study 
of SR and VR in a discrete model with a bichromatic 
input signal and gaussian noise. The model combines the 
features of bistable and threshold setups and turns out 
to be ideal for comparing various aspects of SR and VR. 
Both additive and multiplicative noise are used for driving 
the system in the bistable mode. To our knowledge, this 
is the first instance where VR has been explicitely shown 
in a discrete system. 
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Fig. 16. Variation of response time r for VR and SR as a 
function of noise amplitude E, corresponding to the frequency 
v\ = 0.125, for the cubic map in the threshold setup. For each 
E, the optimum value of the amplitude Z(2) for the higher 
frequency = 0.32 is used. Two regions for VR and SR are 
clearly evident in terms of the response time. 

The system shows bistability in periodic and chaotic 
attractors. Hence a subthreshold signal can be enhanced 
with the inherent chaos in the system in the absence of 
noise. With a composite input signal and noise, all the 
component frequencies are enhanced in the output at dif- 
ferent optimum noise amplitudes. 

Our analysis brings out some fundamental differences 
between the two mechanisms of SR with respect to ampli- 
fication of multisignal inputs. In particular, we find that, 
while the bistable set up responds only to the fundamen- 
tal frequencies present in the input signal, the thresh- 
old mechanism enhances a mixed mode also, a result not 



possible in the context of linear signal processing. This 
may have potential practical applications, especially in the 
study of neuronal mechanisms underlying the detection of 
pitch of complex tones [19,36]. 

Another especially interesting result we have obtained 
is regarding a cross over behavior between VR and SR in 
terms of an initial response time. The response time de- 
creases sharply in going from VR to SR and depends only 
on the background noise level. In the regime of VR, the 
system shows the so called bona fide resonance, where the 
response becomes optimum for a narrow band of interme- 
diate frequencies. Moreover, the optimum value of high 
frequency signal amplitude increases with frequency, but 
as the frequency increases beyond a limit, it becomes less 
effective in producing VR. 

Finally, two other results obtained in the context of 
SR are the usefulness of a high frequency signal in en- 
hancing the SNR of a low frequency signal and the pos- 
sibility of using SR as a filter for the detection or selec- 
tive transmission of fundamental frequencies in a compos- 
ite signal using a bistable nonlinear medium. The former 
is evident from Fig. 8 and shows certain cooperative be- 
havior between the two signals. Similar results have been 
obtained earlier [20,37] under other specific dynamical se- 
tups, but the simplicity of our model suggests that the 
results could be true in general. The latter result arises 
due to the fact that the noise amplitudes for the optimum 
SNR for the two frequencies are different (see Fig. 5 and 
Fig. 8) and the difference AE tend to increase with the 
difference in frequencies (\v\ —V2\). This suggests that SR 
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can, in principle, be used as an effective tool for signal 
detection/transmission in noisy environments. A similar 
idea has been proposed recently [18] in connection with 
the signal propogation along a one dimensional chain of 
coupled over damped oscillators. There it was shown that 
noise can be used to select the harmonic components pro- 
pogated with higher efficiency along the chain. 

The authors thank the hospitality and computing facilities in 
IUCAA, Pune. 
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